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1.0.  INTRODUCTION 


This  report,  prepared  by  the  Analytical  and  Riysical  Simulation  Branch  of 
the  Tank-Automotive  Technology  Directorate,  describes  an  empirical 
approach  to  modeling  the  hydraulic/actuator  portion  of  a  physical  simu¬ 
lation  system.  A  frequency  response  curve  fit  technique  is  applied  to 
measure  actuator  responses.  The  curve  fit  technique  will  not  produce  a 
general  analytical  model  but  will  provide  an  empirical  model  which 
should  incorporate  the  system’s  performance  characteristics  for  many 
types  of  analyses. 


2.0.  OBJECTIVE 

The  objective  of  this  study  is  to  define  a  means  of  predicting  the  vali¬ 
dity  and  fidelity  of  physical  simulation  testing  before  actual  labora¬ 
tory  test  runs  are  made. 

In  the  past,  the  hydraulic/actuator  system  was  assumed  to  be  ideal 
where  the  actuator  responds  perfectly  to  the  desired  command  signal. 
Motions  were  directly  predicted  from  the  command  signals  or  the  source 
from  which  the  command  signals  were  derived,  but  these  prediction  tech¬ 
niques  did  not  account  for  the  performance  capabilities  of  the 
hydraulic /actuator  system.  The  process  used  in  the  past  gave  inac¬ 
curate  predictions  of  accelerations  for  simulation  tests  which  con¬ 
tained  high-frequency  components  (i.e.,  frequencies  beyond  the  system's 
bandwidth) . 


3.0.  CONCLUSIONS 

The  model  technique  discussed  in  this  report  predicts  the  accelerations 
of  the  actuator  with  a  good  degree  of  accuracy.  The  measured  data 
shows  slightly  higher  accelerations  than  the  model  prediction.  A 
Power  Spectral  Density  (PSD)  analysis  illustrates  that  these  differences 
are  caused  by  noise  in  the  measured  data.  The  measured  accelerations 
have  higher  frequency  content  than  the  system  is  capable  of  producing. 


4.0.  RECOMMENDATIONS 

Improvement  in  the  empirical  model  analysis  should  be  made  by  conducting 
a  separate  test  plan  independently  of  the  truck  test  and  by  using 
better  quality  accelerometers. 

The  results  of  this  analysis  could  lead  to  studies  made  on  improving  the 
physical  simulation  performance  beyond  the  current  capabilities.  For 
example,  a  process  could  be  developed  which  would  enhance  the  frequency 
content  of  command  signals  to  compensate  for  the  limited  bandwidth  of 
the  system.  This  would  be  beneficial  for  test  programs  in  which  high- 
frequency  motion  is  essential,  such  as  track  and  suspension  tests.  The 
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truck  test  presented  here  did  not  require  high-frequency  motion  pri¬ 
marily  because  the  truck  cab  was  the  area  of  concern.  Any  high- 
frequency  motion  in  the  truck  wheel  spindles  was  absorbed  by  the 
suspension  system,  resulting  in  negligible  effects  on  truck  cab  motions. 


5.0.  DISCUSSION 
5.1.  General 

The  empirical  model  derived  in  this  report  uses  a  frequency  response 
curve  fit  technique.  It  is  not  considered  a  general  analytical  model  in 
any  sense. 

A  general  analytical  model  would  require  deriving  equations  for  the 
entire  system  by  physical  relationships.  The  curve  fit  method  is  much 
more  easily  obtained  because  it  is  based  on  measured  data  and  does  not 
rely  on  unknown  parameters  which  often  hamper  analytical  modeling.  A 
general  analytical  model  is  more  powerful  for  many  applications  due  to 
the  fact  that  system  alterations  cannot  easily  be  incorporated  into  an 
empirical  model.  The  curve  fit  technique  presented  here  accurately 
simulates  the  response  for  the  linear  operating  regions  of  a  system. 

This  method  is  a  good  starting  point  for  a  general  analytical  model 
development  effort.  The  development  of  the  hydraulic  portion  of  the 
more  general  analytical  model  was  started  during  the  writing  of  this 
report. 

This  work  is  actually  a  part  of  the  physical  simulation  validation 
effort  and  is  independent  of  studies  validating  the  analytical  vehicle 
models  used  to  produce  command  signals  for  laboratory  testing. 

The  results  of  this  study  can  be  used  to  better  describe  the  physical 
simulation  performance  and  can  also  be  used  to  improve  the  prediction  of 
the  actuator  motions  for  a  given  command  signal.  For  example,  actuator 
accelerations  can  be  accurately  predicted  with  a  hydraulic/actuator 
model.  Simulating  the  hydraulic/actuator  system  response  beforehand  can 
save  test  runs  and  may  point  to  potential  problems  with  the  test.  By 
using  a  math  model,  signal  conditioning  requirements  can  also  be  pre¬ 
determined  before  actual  test  runs  are  conducted. 

This  empirical  model  study  was  conducted  during  the  actual  physical 
testing  of  a  5-ton  truck  so  that  actual  measured  test  data  could  be  com¬ 
pared  to  model  results.  Thus,  analysis  and  data  presented  here  are  from 
the  5-ton  truck  laboratory  testing.  The  5-ton  truck  test  was  conducted 
to  evaluate  the  structural  integrity  of  the  truck  cab  and  machine  gun 
mount  subject  to  dynamic  forces  induced  when  the  vehicle  traverses  typi¬ 
cal  cross-country  terrain  profiles.  The  machine  gun  mount  is  known  as 
the  MK19  M0D3  grenade  machine  gun  support  kit  and  has  been  through  many 
design  modifications. 

In  addition,  similar  analysis  was  conducted  with  respect  to  Tank  Test 
Bed  Autoloader  testing.  The  comparison  of  the  model  output  and  tests 
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showed  the  desired  motion  to  be  identical  to  the  actuator  response.  The 
agreement  between  model  prediction  and  actual  tests  is  not  surprising 
because  frequency  content  in  the  turret  region  of  a  heavy  armored  tank 
is  much  below  the  hydraulic  system  bandwidth.  Thus  the  hydrauilc 
actuator  system  easily  responded  to  the  desired  motions  (command 
signals)  during  the  autoloader  tests. 

(Note:  This  report  documents  the  first  attempt  to  model  the  physical 
simulation  tests  and  can  be  used  as  a  reference  for  anyone  interested  in 
doing  this  for  future  test  projects.  It  is  assumed  that  the  reader  has 
some  fundamental  background  in  system  theory  and  is  familiar  with  the 
physical  simulation  techniques  used  in  laboratory  testing  which  are 
described  in  two  reports:  "M9  Driver’s  Hatch  Simulation"!  and 
"Simulation  Test  of  the  MK19  M0D3  Grenade  Machine  Gun  Support  Kit."2) 

5.2.  Procedure 


The  procedure  used  to  derive  a  model  from  frequency  response  data  can  be 
considered  similar  to  the  "black  box  experiment"  conducted  in  a  fun¬ 
damental  systems  course.  The  system  is  treated  like  a  black  box  where 
nothing  is  known  of  the  system  and  the  only  access  to  the  system  is 
through  its  input  and  output  signals.  With  the  assumption  that  the 
system  is  linear,  a  math  model  can  be  derived  by  performing  a  frequency 
response  measurement  of  the  system.  From  the  frequency  response  data,  a 
transfer  function  is  approximated,  which  describes  the  overall  frequency 
characteristics.  Once  the  transfer  function  is  derived,  a  math  model 
can  easily  be  obtained.  This  procedure  was  used  to  model  the 
hydraulic/actuator  system  used  for  the  5-ton  truck  test. 

5.3.  System  Configuration 


The  empirical  model  developed  describes  the  system  configuration 
used  for  the  5-ton  truck  test.  Shown  in  Figure  5-1  is  the  testing  and 
measurement  process  which  was  used  for  both  the  subject  of  the  modeling 
and  to  validate  the  model  results.  A  Computer  Automated  Measurement  and 
Control  (CAMAC)  system  was  used  to  create  the  actuator  position  command 
signals,  and  also  to  sample  the  actuator  accelerations  for  further  anal¬ 
ysis.  This  was  accomplished  using  D/A  and  A/D  converters  respectively. 
The  command  signals  were  created  from  analytical  models  of  the  truck 
traversing  cross-country  terrains.  Their  description  is  beyond  the 
scope  of  this  report.  Further  details  on  how  these  command  signals  were 
generated  are  described  in  the  report  "Simulation  Test  of  the  MK19  M0D3 
Grenade  Machine  Gun  Support  Kit. "3  in  this  case,  the  command  signals 
should  only  be  considered  as  position  command  signals  representing  the 
desired  actuator  motion.  When  conducting  tracked  vehicle  testing,  the 
corresponding  command  signals  are  representative  of  a  terrain  time 
history  at  a  given  speed  and  are  not  dependent  on  analytical  vehicle 
models. 

The  command  signals  generated  are  sent  to  a  signal  conditioner  con¬ 
sisting  of  a  filter  box  which  was  programmed  as  a  5th  order  elliptical 
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5-1.  Testing  and  Measurement  Process 


Produces  simulated 
control  signals 


low-pass  filter  at  6.5  Hz.  This  filter  is  more  dominant  in  attenuating 
high  frequencies  than  the  actuator  system  itself  and  for  our  case  will 
be  considered  as  part  of  the  system  to  be  modeled.  (The  reasons  for 
having  this  filter  in  the  test  are  also  explained  in  the  report 
"Simulation  Test  of  the  MK19  M0D3  Grenade  Machine  Gun  Support  Kit."1*) 

The  filtered  command  signals  are  then  sent  to  the  servoamp  controller 
which  drives  the  actuators.  The  electrohydraulic  system  consists  of  a 
three-stage  servo  valve  system  with  position  feedback  control  for  each 
actuator.  The  position  sensors  used  for  the  main  feedback  loop  are 
Linear  Variable  Differential  Transformers  (LVDTs),  which  measure  the 
displacement  of  the  actuators. 

There  are  six  separate  signals  that  drive  six  independent  actuators  in 
the  vertical  direction  for  the  truck  test.  The  assumption  was  made  that 
the  actuators,  which  are  fastened  to  the  truck  spindle  fixtures,  behave 
independently.  We  are  deriving  a  simple  model,  where  only  one  command 
signal  and  hydraulic  actuator  was  modeled  to  demonstrate  its  capabili¬ 
ties  at  this  time. 

5.4.  Frequency  Response  Curve  Fit 


Both  the  filter  (Signal  Conditioner)  and  the  entire  hydraulic/actuator 
system  will  be  modeled  separately.  Figure  5-2  shows  the  frequency 
response  for  a  typical  hydraulic/actuator  system  used  in  the  physical 
simulation  laboratory.  The  bandwidth  is  in  the  6-Hz  range  (-  3  Db). 

This  response  is  considered  the  closed  loop  response  and  is  obtained  by 
the  input  signal  being  put  into  the  servoamp  controller  and  the  output 
being  the  actuator  position  measured  by  the  LVDT.  The  peak  at  21  Hz  is 
believed  to  be  the  pilot  valve  resonance.  This  measurement  was  made 
using  a  white  noise-random  signal  as  input.  Various  forms  of  frequency 
response  measurements  were  made  which  resulted  in  negligible  differen¬ 
ces.  As  an  example,  several  different  reasonable  levels  of  input  com¬ 
mands  resulted  in  the  same  frequency  response.  This  describes  a  system 
which  is  quite  linear  within  its  operating  range. 

Figure  5-2  also  shows  the  results  of  the  frequency  curve  fit  model. 

There  are  many  frequency  analyzers  and  curve  fit  programs  available. 

The  measurements  and  curve  fit  processes  presented  here  were  done  on  a 
Hewlett  Paokard-3562A  Signal  Analyzer  and  are  described  in  detail  in 
"3562A  Dynamic  Signal  Analyzer -Operating  Manual. "5 

Figure  5-3  shows  the  same  analysis  for  the  filter  system.  The  phase  was 
also  comparable  between  the  measured  frequency  response  and  curve  fit 
model.  The  curve  fit  process  is  based  on  a  weighting  sequence  which  was 
set  to  be  distributed  evenly  for  the  frequency  range  shown.  The  number 
of  poles  and  zeroes  required  for  the  fit  may  be  debatable.  However,  it 
is  apparent  from  the  results  that  the  poles  and  zeroes  in  Table  5-1  are 
more  than  sufficient  for  this  type  of  analysis. 

5.5.  Empirical  Math  Model  Development 


The  poles  and  zeroes  obtained  from  the  frequency  response  curve  fit  are 
now  formulated  into  an  empirical  math  model  for  simulation.  Figure  5-4 
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Table  5-1. 

Frequency  Response  Curve  Fit 

Hydraulic/Actuator  System  (Unloaded) 

Lead  Coefficient 

-50323.04 

Real 

Imaginary 

Zeroes 

226.9 

-4.7 

+/-  144.6 

Poles 

-29.2 

-3.7 

+/-  140.8 

-532.8 

+/-  358.6 

Filter  -  (Used  as  Signal  Conditioner) 
Lead  Coefficient  51.20 


Real  Imaginary 


Zeroes 

57.6 

+/-  59.3 

5.7 

+/-  69.4 

Poles 

-24.1 

+/-  12.8 

-16.8 

+/-  32.8 

-6.2 

+/-  40.6 

Note:  All  poles  and  zeroes  shown  were  obtained  from  the  Hewlett  Packard 
Analyzer  and  multiplied  by  2‘Tr'for  proper  S  domain  formulation. 
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shows  the  simulation  block  diagram  which  describes  the  model  in  terms  of 
Laplace  transforms.  The  transfer  functions  derived  from  this  process 
for  both  the  filter  and  hydraulic  actuator  system  are  given.  Many  forms 
of  computer  software  are  available  to  simulate  this  type  of  model.  The 
one  chosen  for  this  analysis  is  the  Advanced  Continuous  Simulation 
Language  (ACSL).  Details  on  how  to  use  ACSL  are  contained  in  "Advanced 
Continuous  Simulation  Language  (ACSL)  -  Reference  Manual. "6  (See  the 
addendum  for  a  listing  of  the  ACSL  coding  used  to  simulate  the  system 
model.)  FORTRAN  subroutines  were  added  to  process  input  data  and  output 
data  in  Engineering  Research  Division  (ERD)  format  which  is  commonly 
used  by  the  System  Simulation  and  Technology  Division.  ERD  format  has 
been  selected  as  a  common  data  format  used  for  all  physical  simulation 
and  analysis  software.  One  reason  why  it  was  used  in  this  analysis  was 
to  allow  ACSL  to  read  in  the  same  files  which  were  processed  by  CAMAC 
for  command  signals.  In  fact,  the  D/A  converter  is  simulated  in  ACSL  in 
the  same  manner  as  the  (CAMAC)  command  signals  used  in  the  testing. 

5.6.  Model  and  Test  Data  Comparison 


Figure  5-5  shows  a  typical  command  signal  used  for  the  5-ton  truck  test. 
This  particular  one  is  simulating  the  Churchville  1  course  at  10  mph  and 
represents  the  displacement  of  the  front  right  wheel  spindle.  These 
command  signals  are  reproduced  at  100  samples/3ec.  for  both  the  testing 
and  the  model  simulation.  Figure  5-5  also  shows  the  results  of  the  ACSL 
simulation  which  describes  the  actuator  response  to  the  given  command 
signal.  Note  that  the  higher  frequency  waveforms  and  peaks  are  atten¬ 
uated,  which  is  what  would  be  expected.  Strip  chart  recordings  of  the 
test  data  had  demonstrated  about  the  same  wave  form.  The  second  deriva¬ 
tive  of  these  data  was  taken  to  determine  the  spindle  acceleration  and 
is  shown  in  Figure  5-6.  The  second  derivative  was  derived  separately 
from  ACSL  using  a  rise  over  run  (slope)  relationship.  The  measured 
acceleration  from  the  test  is  also  shown  for  comparison.  The  waveforms 
generally  have  the  same  shape  but  the  measured  data  have  higher  acce¬ 
leration  peaks.  The  acceleroraeters/instruraentation  used  for  this 
measurement  have  apparent  noise  problems  associated  with  them  from  pre¬ 
vious  tests.  Complex  high-frequency  components  have  been  recorded  in 
the  past.  Performance  specifications  of  the  accelerometers  are  not 
known  in  regards  to  noise.  Table  5-2  shows  the  results  of  the  study 
in  terms  of  RMS  values  for  the  entire  time  histories. 

The  measured  data  have  higher  magnitudes  than  the  simulated  results  for 
most  cases.  This  is  especially  true  for  the  peak  values.  It  is 
strongly  believed  that  the  noise  of  the  accelerometers  has  created  this 
discrepancy.  For  a  better  comparison,  a  PSD  analysis  was  conducted  and 
the  results  are  illustrated  in  Figures  5-7,  5-8  and  5-9.  These  PSDs 
were  obtained  from  a  computer  program  which  directly  evaluates  the  time 
history  and  transforms  it  into  Fourier  coefficients.  No  window 
averaging  techniques  were  used.  Zeroes  were  added  to  the  end  of  the 
data  so  that  the  number  of  points  was  equal  to  a  power  of  2  for  calcula¬ 
tion  purposes.  This  may  change  the  magnitudes  but  will  not  affect  the 
comparisons  being  made.  The  simulation  and  test  data  were  digitized  at 
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ACTUATOR  POSITION  COMMAND  SIGNAL 
1.50 
1.20 
0.90 
0.60 

IN  0.3O 
0.00 
-0.30 
-0.60 
-0.90 
-1.20 
-1.50 

0.00  2.50  5.00  7.50  10.00  12.50  15.00  17.50  20.00  22.50  25.00 

TIME  (SEC) 


CHU1  10RPH  4  UATTS  UT 
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Table  5-2.  Measured  and 
Right  Spindle 

Simulated  Results 
Acceleration  (Gs) 

for  a  5 -Ton 

Truck  Test  Front 

MODEL  PREDICTION 

MEASURED 

CHV1  10  MPH 

RMS 

.16 

.22 

APG  37  20  MPH 

RMS 

.15 

.16 

CHV1  30  HPH 

RMS 

.18 

.16 

APG 9  8  MPH 

RMS 

.18 

.24 

FORT  KNOX  8  MPH 

RMS 

.19 

.19 
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the  same  rate  of  .004  seconds/sample.  The  PSDs  generally  compare  well 
between  the  measured  data  and  model  results.  There  are  slight  discre¬ 
pancies  but  the  major  differences  show  up  above  the  8-Hz  region  where 
the  model  results  continue  to  drop  in  magnitude  as  would  be  expected  at 
the  higher  frequencies.  However,  the  measured  data  reached  a  plateau^ 
where  they  remain  until  the  anti-aliasing  filter  cut  off  frequency  of  30 
Hz  (not  shown)  attenuates  it.  This  points  more  to  noise  on  the  measured 
data  which  seems  to  be  random  in  nature  as  it  exhibits  white  noise 
characteristics . 
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program  TRUCK_ACT_FI LTER 

This  ACSL  program  was  used  to  simulate  the  Filter  /  Actuator 
system  used  on  the  5  ton  truck  testing.  This  program  is  used 
to  run  time  simulations  to  evaluate  the  transient  response  of 
the  system.  The  transfer  function  coefficients  were  derived 
from  frequency  response  measurements  applying  a  curve  fit 
algor thm  utilizing  a  Hewlett  Packard-signal  analyser. 

Model  -  35S2A 

This  program  first  reads  in  a  ERD  formatted  file  which  is  the 
same  files  used  in  the  testing  (Subroutine  READERD) . 

A  discrete  portion  of  the  program  simulates  the  D/A  converter 
used  in  the  test.  The  output  is  sampled  and  written  out  in 
ERD  format  (Subroutine  WRITEERD)  so  further  analysis  can  be 
conducted. 

The  following  is  the  block  diagram  configuration  in  terms  of 
Laplace  Transforms  which  was  used  in  this  simulation: 


FILTER  MODEL: 


ACTUATOR  MODEL : 
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actuator  model: 


cinterval  cint=0.01 

REAL  INPUT( 7000 ) ,OUTPUT< 7000) 

INTEGER  INDEX1 , INDEX2 ,NSAMP ,NMSAMP  *"USED  FOR  SAMPLING  INDEX" 
"  The  following  arrays  " 

"  are  used  to  describe  the  transfer  function  coefficient" 

ARRAY  FTF3N( 3) ,FTF3D(3) ,FTF4D(3) ,FTF4N(3) 

ARRAY  STF3N(3) ,STF3D<3) ,STF2N<2) ,STF2D<2) 

CONSTANT  STF3N=1 .0,9.32, 20330 . 

CONSTANT  STF3D=1 .0,7. 408,19848 . 5 
CONSTANT  STF2N=1 . 0 , -226 . 93 
CONSTANT  STF2D=1 ,0 ,29.2 
CONSTANT  FTF3N=1 .0 ,-115.16,6827.73 
CONSTANT  FTF3D=1 .0 ,33.7,1357.1 
CONSTANT  FTF4N=1 .0 ,-11 .42,4846.62 
CONSTANT  FTF4D=1 .0,12.4 ,1682 . 74 
INITIAL 

INDEX1=0 

INDEX2=0 

NMSAMP=7000  3>" NUMBER  OF  SAMPLES  OF  SAMPLED  OUTPUT" 

“  kkk'k'k-k'k-k'k-k'k'k'k'k'k-k'k'k'k'k-k'k'kk'k-k'k'k'k'kk'k'kk'k-kk'k'k'k-k-k'k-k-k'k-k-k'k-k-k'kk'k-k  " 

"  -kk-k-k-k-k-k-k-kk'k'k-k  READ  IN  COMMAND  DATA  -k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k" 

11  ■kk-k-k-k'k-k-kk'k'k-k-k-k-k-k-k'k'k'k-k-k-k'k-k-k'k'k-k'kk'kkick-k-k-k-k-k-k-k'k-k-k'k'k-kk-k-k-k-k-k-k  " 

PROCEDURAL (NSAMP , INPUT) 

CALL  READER.D  (  NSAMP ,  INPUT=) 

END  £ "  OF  PROCEDURAL" 

"  -k-k-kmkick-k-k’kmk-k‘k-kick-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k.-k-k1ck'k'k-k'k-k'kmk-k-k‘k'k-k-k-k-k-k-k-k-k  " 

END  %"  OF  INITIAL  " 

DISCRETE  SAMPLE 

"  'k-k-k-kk'k'kk-k-kkk-k'k-k-kk'k-k-k-k'k-k-k'k-k-k'k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k'k'k-k.-k-k-k-k-k-k-k-k-k  " 

" kk-k-k-kk  CREATE  THE  COMMAND  SIGNAL  kkkirk-kkk-kk-kkkk-k-k-k-k-k-kkk " 

"******  BY  SIMULATING  A  D/A  CONVERTER  (100  PO I  NTS/SEC  )***" 

INTERVAL  DTSAMP= .01 
I NDEX1 = I NDEX1+1 

I NF1 = I NPUT ( I NDEX1 )  S"  D/A  CONVERTER  OUTPUT  -  INFi" 

"  ■kk-k-k-k-k-k'kk-k’k-k-k-kk-kk'k-kkk-k-k'k-k'k-k'k’k-k'k-k'k-k'k-k'k'k’k-k’k'k-k'k'k-k’k-k-k'k'k’k-k'k-k-k’k  " 

END  $"OF  DISCRETE  SAMP  1" 
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DERIVATIVE 

" ****************  INPUT  :  INF1  OUTPUT  :  INF5  ********* " 
INF3=51 . 2*1 . 34SE-3*CMPXPL( 1 . 348E-3 ,  .  065 , INF1 ,  0  .  ,  0  . ) 
INF4=TRAN( 2 , 2 , FTF3N , FTF3D , INFS) 

INF5=TRAN( 2 , 2 , FTF4N ,FTF4D , INF4) 

"  ********************************************************** 11 
"****************  ACTUATOR  MODEL  :  ************************ ” 
“ ****************  INPUT  :  INFS  OUTPUT  :  OUTS  ********* " 

INS3=~5G  323 . 04*TRAN( 1,1, STF2N , STF2D , INFS) 

INS4=TRAN< 2,2, STF3N , STF3D , INS3) 

QUTS=2 . 425E-6*CMPXPL( 2 . 425E-S , 2 . 584E-3 , INS4 , 0 . , 0 . ) 
“**********************************************************'* 
END  *"OF  DERIVATIVE" 

DISCRETE  SAMPLE 

•'****  SAMPLE  THE  OUTPUT  SIGNAL  OUTS  ********************** " 
u ****  OUTPUT  IS  DISCRETE  ARRAY  OF  OUTS  ***************** c**'1 
INTERVAL  DTSMP= .004 
STEP= .004 
I NDEX2= I NDEX2+1 
OUTPUT ( INDEX2) =GUTS 
END  *"OF  DISCRETE  SAMP  2" 

DERIVATIVE 

termt ( INDEX2  .GE.  NMSAMP) 

END  *"OF  DERIVATIVE" 

TERMINAL 

“*************  WRITE  OUT  SAMPLED  DATA  IN  ERD  FILE  ******* 
CALL  WRITEERD( =STEP , INDEX2 , OUTPUT) 


END  VOF  TERMINAL" 
END  $" OF  PROGRAM" 


SUBROUTINE  READERD(NSAMP , DATA INPUT) 

* 

*  THIS  SUBROUTINE  USED  TO  READ  IN  ERD  DATA  FILE  AS 

*  INPUT  INTO  ACSL . 

* 

*  CALL  AT  "INITIAL"  PART  OF  ACSL 

*  INCLUDE  AT  "END"  OF  ACSL. 

* 

* 

* 

CHARACTER*80  ERD  TITLE , LONG ,TITLE , DUMMY 80 
CHARACTER*64  ERD^FILE , HDR_FI LE , ERD_FI LE_0 , HDR_FI LE_0 
CHARACTER*32  LONG_NAME( 12), DUMMY 3 2 
CHARACTER*! 2  DUMMY 

CHARACTER*S  SHORT_NAME < 1 2 ) , UN I T_NAME (12), XUN I T , DUMMYS 
CHARACTER* 4  ERD , HDR 
CHARACTER*!  COMMA , REPLY 

REAL*4  SCALE ( 12) ,0FFSET(12) ,DATA( 12 , 30000 ) 
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REAL  RMS(12) ,SMEAN(12) , DATAINPUT( 30000 ) 

INTEGER*2  ERD  UNIT ,  HDRJJNIT ,  I  DATA(  12) 

LOGICAL*4  TIME 
LOGICAL*l  NEWCHAN , RECHAN 
DATA  ERD_UNIT , HDR  JJN IT/10 ,11/ 

ERD  =  ' .ERD' 

HDR  =  ' . HDR' 

kkkkkkkkkkkkkkkkkkkkkkkkkkkickkkkkkirkkkkk'kkkk'k'krk'kk'k’k'k'k-k'k'k'k'k'k'k'k 

■k 

■k  Determine  the  name  of  the  input  data  file 
* 

10  WRITE( 5 , 20 ) 

20  F0RMAT(///,'  This  subroutine  will  send  a  channel  of  a  ' 

+'ERD  file',/,'  to  an  ACSL  simulation',/, 

+'  Enter  file  name  to  send  to  ACSL?') 

READ(5,30)  ERD_FILE 
30  FORMAT (A32) 

CALL  STR*TR I M < HDR  F I LE , ERD  F I LE , LENGTH ) 

HDR  F I LE  <  LENGTH+1 : LENGTH+4 )  =  HDR 
ERD__FI  LE(  LENGTH+1 :LENGTH+4)  =  ERD 

* 

■k  Open  the  data  file,  print  the  header  characteristics,  and  determine  if 

*  this  is  the  correct  data  file 

* 

OPEN ( HDRJJN I T , F I LE=HDR_F I LE , FORM= ' FORMATTED ' , 

+  SHARED , STATUS= ' OLD ' , ERR=21 0 ) 

* 

*  Read  the  header  data 

* 

READ (HDRJJN IT, 60)  DUMMY 
60  FORMAT <A1 2) 

READ  (  HDRJJN  I T ,  70  )  ERD_TITLE 
70  FORMAT (A80) 

READ ( HDRJJN I T , 30 )  NCHAN , COMMA , NSAMP , COMMA , NL I NES , COMMA , NB I N , 

+  COMMA ,NBYTE , COMMA , KEYNUM , COMMA , STEP , COMMA , KEYOPT 

80  FORMAT (6(17, A) , El 3. 6, A, 17) 

READ ( HDRJJN I T , 90 )  SCALE ( 1 ) , ( COMMA , SCALE ( L ) , L=2 , NCHAN ) 

READ ( HDRJJN IT ,90)  OFFSET ( 1 ),( COMMA , OFFSET ( L) ,L=2, NCHAN) 

READ ( HDRJJN I T , 1 0  0 )  ( SHORT_NAME( L) , L=1 .NCHAN) 

READ ( HDRJJN I T , 1 1 0 )  ( LONG_NAME ( L ) ,L=1, NCHAN) 

READ ( HDRJJN IT , 100 )  (UNI T_NAME ( L ) ,L=1 , NCHAN) 

90  FORMAT (13(E13.6,A) ) 

100  FORMAT ( 31 ( A8 ) ) 

110  FORMAT ( 7( A32) ) 

* 

■k  Write  out  the  header  information 
* 

WRITE( 5 , 120 )  ERD_TITLE , NCHAN , NSAMP , STEP 
120  FORMAT (// , '  The  title  for  'this  file  is:",/,"  ',A80,//,'  There  are  * 

+',I2,'  channels  of  data.',/,'  There  are  ',17,'  samples  for  each  da 
+ta  channel.',/,'  The  step  size  is  ',F8.5,'  seconds.',//) 

•k 

■k  See  if  there  are  more  than  16  channels  to  plot 

■k 

•k-k-k-k-k-k-kk'kirk-k-k-k-k-krk-k-kk-kirkk-k-k-k-k-k-k-kkk-k-k-k-k-k-krk-k-k-kk-kirk-k-kk-k-k-k-k-kirk-k-k 
IF  (NCHAN  .GT.  16)  THEN 

TYPE* , 'There  are  more  than  16  channels  to  EDIT.' 

TYPE* ,' Please  FORGET  ABOUT  IT' 

CLOSE  (HDRJJN IT) 

STOP 
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*  *  *  *  *  * 


END  IF 


* 

*  Nrite  out  the  additional  descriptor  lines 

* 

IF  (NUNES  .GT.  0)  THEN 

TYPE*, 'The  following  are  the  optional  descriptor  lines:' 

DO  130  L=1 ,NLINES 

READ( HDR  UN IT, 70)  LONG 
WRITE( 5 ,125)  LONG 
125  FORMAT ('  ',AS0) 

130  CONTINUE 

END  IF 

* 

*  Is  this  the  correct  data  file 

* 

MRITE( 5 , 140 ) 

140  F0RMAT(// , '$1 s  this  the  correct  data  file  to  EVALUATE  (y  or  n)? 

+  ') 

READ (5, 150)  REPLY 
150  FORMAT (A) 

IF  (REPLY  .EQ.  'N'  .OR.  REPLY  .EQ.  'n')  THEN 
CLQSE(HDRJJNIT) 

WRITE( 5,160) 

160  FORMAT(/ , '$Do  you  wish  to  look  at  another  file  (y  or  n)?  ') 

READ ( 5 , 150 )  REPLY 

IF  (REPLY  .EQ.  'N'  .OR.  REPLY  .EQ.  'n')  STOP 
GOTO  10 
ENDIF 

Open  the  data  part  of  the  file 

IF  (KEYNUM  .EQ.  5)  THEN 

OPEN( ERDJJNIT , FI LE=ERD_FI LE , STATUS3' OLD' 

+  , SHARED, FORM3' FORMATTED ' ) 

ELSE 

OPEN( ERDJJNIT , FI LE=ERD_FI LE , STATUS3' OLD' 

+  , SHARED , FORM3' UNFORMATTED' ) 

ENDIF 

CLOSE ( HDR  JJNIT) 

Read  the  data 

J  =  0 
170  J=J+1 

IF  (KEYNUM  .EQ.  5)  THEN 


READ( ERDJJNIT , ISO , ERR=220 , END=230 )  ( DAT A ( I , J) , I =1 .NCHAN) 
ISO  FORMAT ( 19( E13 . 6) ) 

ELSE 

IF  (KEYNUM  .EQ.  0)  THEN 

READC ERDJJNIT, ERR=220 ,END=230)  ( I DATA( I ) , 1=1 , NCHAN) 

DO  190  K=l, NCHAN 

DATA( K , J)  3  FLOAT I ( I  DATA ( K ) ) 

190  CONTINUE 

ELSE 

READ ( ERDJJN I T , ERR=220 , END=230 )  ( DATA ( I , J ;  ,  Is  1 , NCHAN ) 
ENDIF 
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END  IF 


GOTO  170 

★ 

★ 

* 

210  TYPE*, "Error  opening  data  file" 

STOP 

* 

220  TYPE*, "Error  reading  data  in  file" 

230  CLOSE ( ERD_UN IT) 

* 

***********  DATA  IS  READ  IN,  NOW  START  EVALUATION  ******* 
* 

***  Convert  UNscaled  and  UNbiased  data  to  proper  values 
DO  6002  I =1 ,NCHAN 

DO  6003  J=1 ,NSAMP 

DAT  A ( I , J )  =  DATA( I , J)  /  SCALE ( I )  +  OFFSET(I) 
6003  CONTINUE 

SCALE(I )=1 . 

OFFSET ( I ) =0 . 

6002  CONTINUE 

******************************************************** 
******************************************************** 
567 
1018 


+ 


1060 
1050 

1032 

********************************************************** 

* 

*  ASK  WHICH  CHANNEL  YOU  DESIRE  TO  SEND  TO  ACSL 

* 

WRITEC  5 , 5538) 

5538  FGRMATC"  Which  channel  do  you  desire  to  send", 

+  ."  to  ACSL") 

READ (5,*) ICHANACSL 

******************************************************** 

DO  1777  ISMP=1,NSAMP 
DATAINPUTC I SMP) =DATA( I CHANACSL , I SMP) 

1777  CONTINUE 

********************************************************** 

RETURN 

* 

END 


WRITE( 5,1018) 

FORMAT (//, "*****  CHANNELS  *******",//) 

DO  1050  JI =1 ,NCHAN 

IF(JI  .EQ.  6  .OR.  JI  .EQ.  12  .OR.  JI  .EQ.  18  .OR. 
JI  .EQ.  24  .OR.  JI  .EQ.  30  .OR.  JI  .EQ.  36)  THEN 
WRITE( 5 ,1032) 

READ<  5 ,*) 

END  IF 

WRITE( 5 , 1060 ) JI ,LONG_NAME( JI ) , UN I T_NAME ( J I ) 

FORMAT (/ , "  CHANNEL  " , I  2 ,/ , IX , A32 ,/ , IX , A8) 

CONTINUE 
WRITEC  5 ,1032) 

FORMAT (/ , "*********  HIT  RETURN  ******** " ) 

READ( 5 ,*) 


S 
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SUBROUTINE  WRITEERD( STEP ,NSAMP , OUTPUT) 

C  This  subroutine  simply  writes  out  data  arrays 

C  in  ERD  f ormat . 

C 

C  CALL  at  TERMINAL  portion  of  ACSL 

*  C  Include  at  end  of  ACSL  listinq 

C 

CHARACTER*®©  ERD_T I TLE ,  LONG  .TITLE,  DLIMMYSO 
.  CHARACTER*S4  ERD_F I LE , HDR_F I LE , ERD_F I LE_Q , HDR_F I LE_G 

CHARACTER*32  LONG_NAME (16), DUMMY 3 2 
CHARACTER*12  DUMMY 

CHARACTER*®  SHORT_NAME  (IS),  UN  I  T_NAME  (IS),  XUN I T  ,  DUMMY'S 
CHAR.ACTER*4  ERD ,  HDR 
CHARACTER*1  COMMA , REPLY 
CHARACTER*2  I  OPERATE (20) 

DIMENSION  SMAX(IQ) ,SMIN(18) 

REAL*4  SCALE(IS) ,0FFSET(16) , DATA( 16 , 30000 ) 

REAL  RMS ( 1 6 ) , SMEAN ( 1 S ) , OUTPUT ( 30 0 0 0 ) 

INTEGER*4  START_.SAMP_.ELIM ,  END_SAMP  ELIM 
INTEGER*2  ERDJJNIT , HDRJJNIT , I DATA( IS) 

LOGI CAL*4  TIME 

LOG I CAL*1  NEWCHAN , RECHAN 

DATA  ERDJJNIT, HDRJJNIT/10, 11/ 

„  ERD  =  ' .ERD' 

HDR  =  ' .HDR'  ■ 

* 

WRITE( 5 ,1831 ) 

1881  FORMAT(///,'  ENTER  how  many  channels  are  output') 

READ ( 5 , * ) NCHAN 

* 

DO  1998  J=l, NCHAN 
WRITE( 5 , 1888) J 

1888  FORMAT (/// , '  ENTER  the  LONG  NAME  for  channel  ',12) 

READ ( 5 , 1889) LQNG_NAME( J) 

1339  FORMAT (A32) 

WRITE( 5 , 1992) J 

1992  FORMAT (//,'  ENTER  the  SHORT  NAME  for  channel  ',12) 

READ <5,1 920 )SHORT_NAME(J) 

1920  FORMAT (AS) 

WRITE( 5 , 1993) J 

1993  FORMAT (// , '  ENTER  the  UNIT  NAME  for  channel  ',12) 

READ ( 5 , 1920 ) UNIT_NAME( J) 

* 

OFFSET ( J) =0 . 

SCALE( J) =1 . 

* 

1998  CONTINUE 
* 

DO  176S  I SMP=1 .NSAMP  ,  ' 

DATA (1,1 SMP) =OUTPUT ( I SMP ) 

176S  CONTINUE 
* 

WRITE( 5,201) 

201  FORMAT (//,'  Indicate  how  new  data  file  is  to  be  stored. 
+'  0=2  byte  integer  (binary)',/, 

+'  1  =  4  byte  floating  point  (binary)',/. 
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+  '  2  =  8  byte  floating  point  ( bi nary ) ' ,/ , 

+  '  3=8  byte  complex  (binary)',/, 

+  '  4  =  IS  byte  complex  (binary)',/, 

+'  5  =  formatted  floating  point.  The  format  is  (Nchannels) E13 . 

+S' ,/, '$Enter  selection  (0-5) :(1  CHOSEN  MOST  COMMONLY)  ') 

READ  (5,*)  KEYNUM 

★ 

*  do  not  let  the  user  choose  complex  numbers 

* 

IF  (KEYNUM  .EQ.  3  .OR.  KEYNUM  .EQ.  4)  THEN 
TYPE*,'  ' 

TYPE* ,' Choose  another  format  besides  complex  numbers.' 

ELSE  I F( KEYNUM  .LT.  0  .OR.  KEYNUM  . GT .  5)  THEN 
TYPE*,'  ' 

TYPE* j 'Selection  out  of  range.' 

END  IF 

* 

*  open  and  create  files 

* 

* 

*  begin  writing  header  information 

* 

HRITE( 5 , 265) 

265  FORMAT (//, '$En ter  name  of  the  data  file  to  write  to:  ', 

+/,'  ERD  FORMAT  ASSUMED') 

READ (5, 267)  ERD_FILE_0 
267  FORMAT (A32) 

*  Create  the  two  file  names 

* 

UIRTTFf®;  4446’) 

4446  FORMAT ('  Enter  ERD  title?') 

READ ( 5 , 4447) ERD_TITLE 

4447  FORMAT (A30) 

HDR_FILE_0  =  ERD_FILE_0 

CALL  STR*TRIM(HDR_FILE_0,ERD  FI LE_0 , LENGTH) 

HDR  FILE_0(LENGTH+1 :LENGTH+4T  =  HDR 
ERD_FI LE_0( LENGTH+1 :LENGTH+4)  =  ERD 
0PEN(  HDR_UNIT  ,  FI  LE=HDR  FI  LE__0 ,  STATUS*'  UNKNOWN'  , 

+  FORM*' FORMATTED' ,RECL=25S) 

* 

*  WRITE  OUT  HEADER  DATA 

* 

* 

*  UNKNOWN  KNOWNS 

* 

DUMMY  =  ' ERDF I  LEVI .00' 

KEYOPT =0 

NLINES=0 

NBIN=-1 

NBIN=-1 

NBYTE=-1 

COMMA*' , ' 

WRITE( HDR_UNIT , 270 ) -DUMMY 
270  FORMAT ( A12) 

WR I TE ( HDR  JJN IT , 280 )  ERD_TITLE 
280  FORMAT (A80) 

WR I TE ( HDR JJN I T , 280 )  NCHAN , COMMA , NSAMP , COMMA , NL I NES , COMMA , NE I N , 

&  C  OMM A  ,  N  BYT  E  ,  C  OMMA  .  K  EYN  UM  ,  C  OMM A  ,  ST  E  P  ,  C  OMMA' ,  K  E Y  0  FT 

230  FORMAT ( 6 ( I  7 , A ) , El 3 . 6 , A , I  7 ) 

WR I TE ( HDR  JJN I T , 30  0 )  SCALE (!),(( COMMA , SCALE ( J) ) , J=2 , NCHAN ) 
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300  FORMAT (18(E13.6,A)) 

WR I TE ( HDR_UN I T , 30 0 )  OFFSET(l) ,(( COMMA , OFFSET( J) ) ,J=2,NCHAN) 
WRITE(HDR_UNIT,310)  ( SHORT_NAME( J) ,J=1,NCHAN) 

310  FORMAT ( 31 (  A8) ) 

WRITE( HDRJJNIT ,320)  ( LONG_NAME( J) ,J=1,NCHAN) 

320  FORMAT (7(A32) ) 

WR I  TE  ( HDR_UN I T , 31 0 )  (  UN  I TJMAME  (  J  > ,J=1,NCHAN) 

* 

*  write  9+  lines  to  file 

* 

C  TYPE*,'  ' 

I F(NLINES  .EQ.  0)  GOTO  330 

*  TYPE*, 'Enter  additional  descriptor  lines' 

DO  330  J=1 ,NLINES 

*  READ (5, 280)  LONG 

WRITE(HDR  UNIT, 6560) J , SHORT_NAME( J) , RMS( J) ,SMEAN( J) 

6560  FORMAT ( '  CHAN  ' , 1 2 , 3X , A3 , 3X , '  RM3=  ',E15.3,4X,'  MEAN=  ',E15.3) 
330  CONTINUE 

* 

CLOSE (HDRJJN IT) 

*  write  data  to  file 
■k 

IF  (KEYNUM  .EQ.  5)  THEN 

OPEN ( ERD_UN I T , F I LE=ERD_F I LE_0 , FORM= ' FORMATTED ' , 

+  STATUS= ' UNKNOWN ' ,  RECL=256) 

ELSE 

OPEN ( ERD  JJN I T , F I LE=ERD_F I LE_0 , FORM= ' UNFORMATTED ' , 

+  STATUS=' UNKNOWN' ,  RECL=256) 

END  IF 

*  2  byte  integer 

* 

IF  (KEYNUM  .EQ.  0)  THEN 
DO  350  J=1 ,NSAMP 
WRITE(ERD_UNIT,ERR=406) 

+  ( 1 1 FIX( DATA( L , J) )  , L=1 ,NCHAN) 

350  CONTINUE 

•k 

*  4  byte  floation  -  binary 

i k 

ELSE  IF  (KEYNUM  .EQ.  1)  THEN 
DO  370  J=1 ,NSAMP 

WRITE(ERD_UNIT,ERR=406) 

+  (DATA(L, J) ,L=1 ,NCHAN) 

370  CONTINUE 

★ 

*  3  byte  floating  -  binary 
■k 

ELSE  I F( KEYNUM  .EQ.  2)  THEN 
DO  330  J=1,NSAMP 
WRITE( ERD  JJN IT , ERR=406) 

+  (DBLE(DATA(L, J) ) ,L=1 ,NCHAN) 

390  CONTINUE 

■k 

*  formatted  output 
•k 

ELSE 

DO  410  J=1 ,NSAMP 

WR I TE ( ERD  JJN I T , 40  5 , ERR=40  6 ) 

+  ( DATA ( L , J) , L=1 , NCHAN ) 
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410  CONTINUE 

405  FORMAT ( 19( E13 . 6) ) 
GOTO  35 

406  WRITEC  5 ,407)  J-l 

407  FORMAT </,'  There 
+le  filled  up . ' ,/, '  C 

ENDIF 

35  CLOSE (ERDJJNIT) 

■Jr 

RETURN 

★ 

END 


,407)  J-l  ,  .  .  . 

/,'  There  were  ',110,'  records  written  out  before  the  ri 

up.',/,'  Change  NSAMP  in  the  header  file  accordingly . ' ) 


* 


* 
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